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Q ■ Abstract 

The Bousso-Polchinski (BP) Landscape is a proposal for solving the Cosmological 
Constant Problem. The solution requires counting the states in a very thin shell 
^ ■ in flux space. We find an exact formula for this counting problem which has two 

simple asymptotic regimes, one of them being the method of counting low A states 
_ given originally by Bousso and Polchinski. We finally give some applications of the 

\^ • extended formula: a robust property of the Landscape which can be identified with 



an effective occupation number, an estimator for the minimum cosmological constant 
' and a possible influence on the KKLT stabilization mechanism. 
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1 Introduction 

The eternal inflation picture of the multiverse consists of de Sitter bubbles nucleating in 
certain vacuum state of very high energy density [H El El IH E]. Bubbles can be created 
inside other bubbles, and this provides a dynamical relaxation mechanism which gives rise 
to an average neutralization of the cosmological constant [6l [7]. We may wonder if it's 
possible to formulate a model of eternal inflation with relaxation containing vacuum states 
of a cosmological constant as small as the observed valu^ [HI [9] 

Aob. = 1.5 X 10-^2^ (1) 

without fine-tuning the parameters of the model and avoiding the need of invoking the 
anthropic principle [101 [HI [12] , which has been used to explain why we are living in such 
a special region of the multiverse. The smallness of the number in ([T]) is the cosmological 
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constant problem [131 [13] • An attempt for a solution is given in the Bousso-Polchinski 
Landscape [15], in which a large amount J of quantized fluxes of charges ... j leads 

to an effective cosmological constant 

1 ^ 

A = Ao+2E^M- (2) 

In ([2]), Aq is a negative number of order —1, and the integer J-tuple (ni, ■ ■ ■ ,nj) charac- 
terizes each of the vacua of the Landscape. Without fine-tuning, for large J and incom- 
mensurate charges {qj} this model contains states of small A. The problem arises now as 
how to count them. Unfortunately, the amount of these anthropic states is expected to be 
very small as compared to the total number of vacua in the Landscape [17] , and therefore 
it seems not to be other way out but invoking the anthropic principle to explain the value 

The states in the Bousso-Polchinski Landscape can be viewed as nodes of a lattice in 
flux space R'^. We call this lattice C and the charges qi are the periods of C, that is, 

£ = {{riiqi, ■ ■ ■ ,njqj) E R'' : rii, ■ ■ ■ ,nj E Z] . (3) 

A single state A in this lattice is characterized by the quantum numbers (ni, ■ ■ ■ , rij), and 
its cosmological constant is, according (|2]), 

A(A) = Ao + ^||Af . (4) 

Vacuum states in the Bousso-Polchinski Landscape are defined in the semiclassical approx- 
imation as stationary points of an effective action. If we consider two neighbor states of 
very high A (neighbor states have quantum numbers {fij} which are different at one place 
and by one unit) we find that the energy barrier separating them is small, that is, the me- 
diating Brown- Teitelboim instanton has a comparatively low action. These states are not 
isolated and consequently the semiclassical approximation breaks down for them. More- 
over, when A reaches a value near the Planck energy density (of order A 1), quantum 
gravity effects become important, and some approximations made in the Bousso-Polchinski 
model (as neglecting the backreaction effect, for example) are no longer valid. 

Thus the Bousso-Polchinski Landscape is a finite subset (yet an enormous one, a com- 
monly quoted number being 10^°° [HI [HI [201 [21]) of the lattice (jSj), comprising the nodes 
with cosmological constant smaller than some value Ai = 0{1). 

We will review the counting argument of Bousso and Polchinski [151 [iH]- Around each 
node A of the lattice C we have a Voronoi cell which is a translate of the parallelotope Q = 
Y[i=i[~^^ f"] of volume vol Q = Y[i=i Qi- On the other hand, each value of the cosmological 
constant Aq < A < Ai defines a ball B'^{Ra) in flux space of radius -Ra = a/2(A — Aq), 
and whose volume is ^ 

YolB-'iRA) = ^YO\S^~\ (5) 
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where the volume of the J — 1 sphere is 



The BP counting argument consists of computing the number of states inside a ball of any 
radius r, which will be called Qj{r), by taking the quotient between the volume of the ball 
and the volume of the cell of a single state: 

vol B'^ (r) 

Plugging in some numbers, if Aq = —1, Ai = 1, J = 300 and qi = ^, we have a BP 
Landscape with 

(^300(2) = 1.3 X 102°2 states. (8) 

This is certainly very small when compared with 10^°°. But if we take a charge ten times 
smaller = 0.01, we obtain f23oo(2) = 1.3 x 10^°^. Nevertheless, there is another argument 
which gives a more impressive number. The quotient between 2i?Ai and qi gives the number 
of nodes that fit in axis i. So, the total number of nodes is the product of all these quotients, 
which is also the volume of the hypercube of side 2R\^ divided by the volume of the cell. 
With the same numbers this quantity is 4.1 x 10^*°. This argument is wrong though, 
because the vast majority of the nodes of the hypercube lie outside the sphere of radius 
i?Ai, and thus they are not states of the Landscape. 

We may wonder how many states there are with values of A comprised between and a 
small positive value of the cosmological constant A^. Following the same argument as in the 
preceding paragraph, we compute the quotient between the volume of the shell comprised 
between radii Ro = a/2|Ao| and = a/2(A7^^A^ ^ Rq + ^ and we obtain 

vol S''^~i 

K = ^ARe) - ^ARo) ^ Rt' — VTT • (9) 

volQ 

If Ae is of the order of the observed value ([T]), with the previous numbers we obtain 
A/'e = 2.1 X 10^^. Thus, if the degeneracy of these states is smaller than this number, we 
can find states with a realistic cosmological constant. 

Nevertheless, as the authors of [15| point out, this argument is not valid when any of 
the charges qi exceed Rq/\/~J. We can see that strange things happen as J grows with all 
the charges fixed. There is a critical value of J above which the volume of the sphere is 
smaller than the volume of a single cell, 

vol;B'^(/2o) < volQ, (10) 

and thus computing their quotient is not useful for counting. From inequality (flQj) we have 

^ > 27re > 17 (11) 

Rq 
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with q = -(/vol Q. Other strange thing happen for large J. Let's assume for simphcity 
that all charges are equal to q. The corner of the cell centered at the origin is located at a 
distance f from it, and it reaches (and surpasses) the radius Rq when 



Jg2 



>4, 



(12) 




so we can expect the angular region near the corner to be devoid of states. Thus, the 
distribution of low A states is not isotropic in flux space, and in particular, it has no 
spherical symmetry. 

All these conditions coincide: when the parameter h = ^ is large we cannot count by 
dividing volume^. But there are instances of the BP Landscape in which the parameter h 
can be large and nevertheless the model contains a huge amount of states. In such cases the 
formula ([7]) should not be used, and another formula is needed. Other counting methods 
in the BP Landscape have been proposed so far [161 US l22l EH [24], but all of them have a 
limited range of validity. 

The remainder of the paper is organized as follows. In section [2] we will propose an exact 
counting formula which is reduced to ([7]) for small h. We will also provide an asymptotic 
formula for the regime of large h. In section [3] we extend the method used previously to 
the study of other properties of the BP Landscape, in particular the counting of low-lying 
states, an estimate of the minimum value of the cosmological constant and the possible 
influence of the non-trivial fraction of nonvanishing fluxes in the KKLT moduli stabilization 
mechanism. The conclusions are summarized in section HI 

2 The BP Landscape degeneracy 

In this section, we will obtain an exact integral representation for the number of nodes 
of the lattice inside a sphere of arbitrary radius, and we will analyze its main asymptotic 
regimes. 

2.1 The exact representation 

We start with the number of nodes in the lattice inside a sphere in flux space of radius r. 
This magnitude is called fij(r) above: 



In the previous equation, vertical bars denote cardinality. An alternative expression can 
be given in terms of the characteristic function of an interval / 



nj(r) = \{\ e C: \\\\\ < r 



} 



(13) 




(14) 



^With the numbers given above, h — 0.75, so we can trust eq. (O. 
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so that 

njir) = J2x[o,r]m\)- (15) 

Expression ( |T5i) is exact, and the sum is extended to the full lattice, whitout any problem 
given that x function adds 1 for each node inside the sphere, and therefore the result is 
always finite. Clearly, ( |T5i) is equivalent to directly counting the nodes (the "brute-force" 
counting method), hence it cannot be used in order to obtain numbers as in ([7]). 
The density of states associated to (fT5|) is 

^j{r) = ■ (16) 

which will be called the "BP Landscape degeneracy" . By writing the characteristic function 
in terms of the Heaviside step function 

X[o,.](||A||) = e(||A||)-^^(||Af -r^), (17) 

we obtain 

ujj{r) = 2rY,S{r^ -Uf) ■ (18) 

The counting function ^lj{r) is a stepwise monotonically non-decreasing function, and thus 
its derivative u}j{r) is a sum of Dirac deltas. It is supported at those values of r which 
correspond to the values that are actually attained by the norms of the lattice nodes. Let 
Ai be the set of these values; we have 

Uj{r) = 2rY,^jif^W-f^') (19) 

where we have defined the true degeneracy wj{fi) as the integer- valued function which 
counts the number of nodes in the lattice C whose norm is /i, that is, the number of 
decompositions of a number as a sum X]j=i 1]^% for ni, ■ ■ ■ ,nj integers and qi, - ■ ■ ,qj 
arbitrary real numbers. 

We can express the Dirac delta which appears in flTSl) as a contour integral: 

5p_||Af) = J- / e^(^'-ll^ll')ds, (20) 

where the contour is a vertical line crossing the positive real axis, 

J = {c + iT : T eR,c> 0} . (21) 
Substituting ([20]) in ([H]), we obtain 

2r 



.xec 



e-llA| 



ds . (22) 
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This particular representation allows us to perform the sum extended to the whole lattice: 

J 



Mr) = I- / e" 



2r 
2r 



E En 

.711 "jezi=i 



ds 



■ J 



ds 



ds. 



The sum is hidden in the function 



(23) 



(24) 



valid for Re s > 0, which is a particular case of a Jacobi theta function: 

^3(^;g) = 5^g"'e2™% 



(25) 



for complex z and q with |g| < It satisfies the functional equation 

^2, 



n& ' meZ * 



which is a consequence of the Poisson summation formula. 

Our exact formula for the BP Landscape degeneracy is then 



which is an inverse Laplace transform, that is, 



ds - 



(26) 



(27) 



e uj{r\ 



dr = JJ^9(sg|). 

3=1 



The integration of (1271) with the initial condition fij(O) = 1 gives 



f^j(r) = 1 + 



1 

2Td 
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ds, 



(2J 



(29) 



■^The second argument q of Jacobi theta functions, the so-called nome, shouldn't be mistaken with the 
charge q. 
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We will close this subsection with a final remark. By Laplace transforming (fT9!) and 
comparing it with its alternative form fl28|) we obtain 

J 

Y^WJ{^,)e-^^' = \{d{sq^). (30) 

Substituting all qj = 1, we can see that the possible values of the numbers /i^ when jj, E M. 
are those which can be represented as the sum of J integer squares. These are all the 
non-negative integers, obtaining in this case 

oo 

l + 2^wj{n)e-''' = d{sy . (31) 

n=0 

Formula f l3T|) is the generating function of the number of different decompositions of a pos- 
itive integer n as the sum of J integer square^. Thus, fl30|) can be taken as a generalization 
of dS]). 



2.2 The large distance (or BP) regime 

Now we will turn to the approximate evaluations of integrals (!27|) and (!29|) . For this 
purpose we need the asymptotic behavior of d function. 

Function d{s) has two simple asymptotic regimes for real and positive s, as can be seen 
from the functional equation ( |26l) : 

^{s) ^ ^ and ^9(s) 1 + 26'' . (32) 

We can visualize these asymptotes by plotting the logarithm of the quotient between and 
each of them. This is done in figure [H where we can see that the limit s — )■ is (reasonably) 
valid for s < 1 and the limit s — )■ oo is valid for s > 2. In the middle regime s G [1,2] none 
of the two former cases is accurate enough, and we will have a mixed, interpolating regime 
between them. 

The first case we will consider is s — )■ 0. In this regime, we simply make the integration 
contour 7 pass near the origin in the complex plane, where has a singularity. Assuming 
that the main contribution to the integral will come from this region, we can replace i} by 
its asymptotic value when s — )■ and write 

2r f 

'^'^''^ " I;' 

■^In number theory, the number of decompositions of a positive integer n as the sum of J squares is 
caUed rj(n) and its generating function is usually written using a variable x = e^* with < 1. 




ds 



(33) 
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Complementary asymptotic approximations of *(s) 




s (real) 



Figure 1: Complementary asymptotic regimes of the i9 function. Logarithms of the quotient 
of I? and its asymptotes are shown. A good approximation can be seen as a flat line at 
zero height. Both regimes cross near 1.5, which is the center of the interval [1,2] where 
the accuracy is lower. Note that the horizontal line at 0.01 signals that the quotient is 1% 
different from 1. All plots in this paper were done using R |25j . 
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This integral is an elementary inverse Laplace transform: 



vra 2r f „2 ds 2tc2 r'^ 



Equation fl34l) is the derivative of ([7]), that is, BP count. It is valid for large r distances, 
because it has been obtained using small s region (note that a well known property of 
Laplace transform pairs is that the asymptotic behaviour of a signal for large r is deter- 
mined by the small s behaviour of its transform and vice versa). For this reason we call 
the formula the large distance regime, or BP regime. 

But we can make the restriction imposed on the distance more quantitative in the 
following way. The best choice for the real part of the contour is the saddle point of 
integral fl34|) . which is the stationary point of the function 

= sr^ — — logs. (35) 



(36) 



This saddle point is 

. J 
s* = — 

2r 

and in its vicinity we find the most important contribution to the integral. But the re- 
placement of the ^9(s) functions by its small-s behavior is valid only if the argument s is 
less than 1 (see figure [I]), so we must have 



s*q^ = ^<l for all J. (37) 



Nevertheless, condition ( l37ll does not guarantee the validity of the replacement ■}){s) — >■ -^y^ 
on the integration contour away from the real axis. A stronger restriction is imposed by 
demanding the applicability of the saddle point method in this regime. If the steepest 
descent approximation is valid, then the main contribution to the integral comes from the 
vicinity of the saddle point, which justifies the replacement of the asymptote. The exact 
and approximate evaluations of flM|) will have the same validity if both methods reach the 
same result. But the saddle point approximation on f l34l) has the effect of using Stirling's 
approximation on the gamma function. Thus, both methods agree if J is large enough. 

The correctness of the saddle point approximation of ( IMI) can be assessed by rewriting 
it in the form 



ujj{r) ^ vr 2 



f e^'^'-i i°s(5'«) ds with log g = 4 V log q^ • (3^ 
The change of variable q^s = w transforms the exponent of the integrand into 



\'W 1 
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where h = The vahdity condition of the saddle point approximation is (j){w*) ^ 1 
with w* = h/2 the stationary point of 0. This condition is fulfilled if J is large and 

— --logw* = -(l-log-j >1 /i<-^ 0.736, (40) 

which is analogous to the condition stated by Bousso and Polchinski for the validity of 
their formula. We have derived it as a validity condition for the small-s asymptotic regime 
of the exact counting formula^. 

Finally, the large J condition controls the validity of Stirling's approximation for the 
gamma function. But this restriction is not needed because the integral has been done in 
closed form. Thus, only condition ( HOl) remains. 



2.3 The small distance regime 

In this case we are in the regime in which the asymptotic expansion of for large values 
of its argument is valid. We can write flj{r) in f l2^ as 

njir) = l + ^ [ fis)e^^-'^ds, with , , (41) 



27rz 
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(s) = sr2 + ^/^^ log(l + 2e-'?.^) . 



The saddle-point approximation to this integral is given by 

0,M.l + _L.;„.,e«.-.^. (42) 

where s* is the stationary point of 0(s). The saddle point s* is a minimum for real s; 
hence, the steepest descent contour crosses vertically the real axis and coincides locally with 
7. Unfortunately, we cannot solve the saddle point equation in closed form for arbitrary 
charges. Nevertheless, in the simplest case in which all charges are equal qi = ■ ■ ■ = qj = q, 
we obtain 

,2 



r 



^s-q^ +2 ^ 



.V = log2(^-l). (43) 



The saddle point computed through (H3l) is consistent with the regime of large argument 
of ^9 if 

T 2 J 2 2 

s*q^ = \og2{^^-l)> 2 l|_ > 1 + |. ~ 4.694 . (44) 

This condition is satisfied for fixed charge q and dimension J if the distance is small enough; 
for this reason this regime is called the small distance regime. 



^Incidentally, the adimensional parameter h — occurring in (p9| resembles the t'Hooft coupling in 
the so-called planar limit of field theory, in which the number N characterizing the gauge group tends to 



infinity and the Yang-Mills coupling constant ^ym vanishes with the product Ngy^^ (the t'Hooft coupling) 



held fixed. 
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In terms of the parameter h = the approximate saddle point (which will be called 
u{h) = q'^s*{h)) is, gluing together eqs. ( 1371 143|) 



2 _ )lh if /i< 2, 



Eq. fH5l) is plotted in figure [2], along with the numerical solution obtained in a range of h 
which is not covered so far, but will be considered in the next subsection. 




Figure 2: Numerical saddle point u{h) of the ujj{r) integrand for equal charges. For h < 2 
it agrees with its low-h asymptote, and for h > 5 it agrees with its high-/i asymptote. The 
mixed regime, where the numerical solution smoothly interpolates between the asymptotes, 
is shown in the central rectangle. 



Substituting the large h regime for u{h) leads to 

-2 -log2(/i-l) 



(Pis*) = ^ u{h) + Jlog(l + 2e-"('*)) = J 

ill / *\ ^dq c T 4 ' 



[1 + 2e-«W)2 

Ju(h) 

2 1 — e h. 



h 



rlog2(fe-l) 



log2(/i-l) ' 



h 



log 



h-lJi 



(46) 
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and the saddle-point approximations for Qj{r) and wj(r) result in 



1 {2h^2)i -1 f h y+h 
'^'^~^V2^ log(2/.-2) U-J 



(47a) 



ujj{r) = ^ ^^(l 7) • (47b) 



We must note that these magnitudes depend on r through h = It should be stressed 
that f l47b|) has been obtained by using the saddle point approximation of f l7r|) and not by 



differentiating f l47ap : although these two results are asymptotically equivalent, they differ 
in subleading terms. 

We will now consider the validity of formulae f l47a[ I47bp . The saddle-point approxima- 



tion is good if 0(s*) ^ 1. We could rewrite 0(s*) as follows: 

log2(/i-l) / h 



jm with m= ^^^^ 



and, for large J, demand ip{h) > 1 as we did before. But this would contradict f l44|) . so we 
better choose some Jq > 2 and write 

0(s*) = ^Jom , (49) 

-'0 

which satisfies (/)(s*) 1 if J/Jq is large and Joip(h) > 1. For example, we have 

Jo = {1,2,3,4,5} /i< {2.2,7.2, 12.5,18.3,24.3,30.6}. (50) 

This restriction implies that the distance cannot be too small in order to preserve the 
validity of the saddle-point approximation. 

For very small r, we can choose a large real part of s and approximate (1 + 2e^^ '^)-^ ^ 
1 + 2Je~'^ and then flj{r) is reduced to 

nj{r) ^e{r^) + 2je{r^ -q^), (51) 

that is, only the node at the origin and its 2J neighbors contribute to Qj. The formula 
fl47al) cannot reproduce this result, and therefore there must exist a validity condition 
which forbids too small distances. This is exemplified in fISUl) . However, that restriction 
turns out to be of no importance because close-to-the-origin nodes represent a negligible 
fraction of the whole. 

2.4 The middle distance regime 

When h takes a value in which no accurate asymptotic approximation of 'd is available, 
we are in the middle distance or crossing regime. In this situation the saddle point can 



12 



be computed only numerically. This has been done in figure [2] by solving the following 
equation 



whose solution u{h) is the saddle point. This solution can always be obtained, with no 
condition on the value of h, but it coincides with (l45ll in the specified ranges. Using this 
solution, the saddle-point approximation of ujj{r) can be computed, but only for large J, 
and only for not too small distances. 

In figure 131 we compare equations flMj) and fl47b|) . displaying between them the crossing 
regime. Both small and large distance regimes show up for constant J and q at different 
distances r. Note that small h corresponds to large r and vice versa, but "small" and 
"large" distances are J-dependent concepts, so that both regimes can have their own range 
of validity. Moreover, for sufficiently high J almost all relevant distances in flux space can 
be considered "small". In such cases, the BP formula (17fM|) should be replaced by the 
correct asymptotic one given in fl47af47bp . 

In the case of different charges in the small distance regime, the approximate saddle- 
point equation cannot be exactly solved. Thus, one needs to solve the complete saddle-point 
equation 



The h parameter does not appear in this equation, and it is not clear what kind of average 
charge must be used to define it. 

However, there is a class of models in which the crossing regime dominates over the 
small and large distance regimes. It is enough to consider a charge distribution in which 
the smallest and biggest charges are well separated. In such cases, there will be regions 
in the s plane where every ■(9 factor could in principle lie in a different asymptotic regime, 
and thus we will obtain a plethora of intermediate regimes in which neither BP nor small 
distance regimes will be accurate enough. For large J, the degeneracy density uj{r) can 
be obtained by computing the numerical solution of ([S2D and then using it in the saddle 
point approximation of the exact integral (|27|) . 

3 Applications 

In this section, we will show how uj{r) helps to estimate other properties of BP models. 
We will consider the number of states in the anthropic window, the distribution of non- 
vanishing fluxes of a typical state, the minimum cosmological constant and a possible 
consequence on the KKLT moduli stabilization mechanism. 




(53) 
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Degeneracy (J = 10, q = 0.15) 



Degeneracy (J = 150, q = 0.15) 




Degeneracy (J = 300, q = 0.15) Degeneracy (J = 900, q = 0.15) 




n I I I r n r i i r 

0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 



Figure 3: Asymptotic regimes of the Landscape degeneracy ooj{r) for equal charges q = 0.15 
and different J as functions of r. The actual ooj{r) is close to the upper envelope of both 
curves (green dashed line). For small J the large distance or BP regime dominates (red line), 
but for large J the small distance regime (blue line) spans the whole [0, 2] interval. Vertical 
dashed lines delimit the crossing regime, corresponding to /i G [2, 5], which is above r = 2 in 
the J = 900 panel. 
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3.1 Number of states in the Weinberg Window 



The number of states of positive cosmological constant bounded by a small value is the 
number of nodes of the lattice in flux space whose distance to the origin lies in the interval 
[Ro, Re], where Rq = a/2|Ao| and Re = \/2{Ke - Aq) ^ -Rq + ^ so that the width of the 
shell is e = 4^: 

Ho 

Me = QjiRe) - ^j{Ro) ~ UJj{R^)^ . (54) 

We remind the reader that fij(r) is the number of states inside a sphere of radius r in flux 
space and uj{r) = If A^ is the width of the anthropic range Aww (the so-called 

Weinberg Window), then the number of states in it is 

(55) 

Kq 

Computation of uj{Ro) should be done along the lines of the previous section. Thus, 
the expression fl55|) can be used for all values of h using the relevant approximation of the 
exact formula fl271) . which includes the BP regime as well as the small distance and crossing 
regimes. 



3.2 Typical number of non-vanishing fluxes 

As before, consider the nodes of the lattice C having cosmological constant between and 
Ae- They will lie inside a thin shell of width e = Re — R above radius R (former radius Rq 
will be called R in this section). The set of nodes inside the shell will be called Eg: 

= {A G £: i? < ||A|| < i?J \j:e\=K- (56) 

We will assume that e is smaller than the charges qi so that (!54l) is valid but A/'e ^ 1. 

Taking J = 2, we will find at most four nodes in the shell with one vanishing component. 
Thus, the remaining states will have two nonzero components. In the J = 3 case, the states 
in the shell are located at the axes (at most six) with only one nonzero component, at the 
coordinate planes with two nonzero components (a larger charge-dependent number) and 
at the "bulk" of the sphere with all three nonzero components (the most abundant). In 
this way, we find that the typical number of non-vanishing components is J for the cases 
J = 2,3. 

Thus, after drawing a node of the shell at random (assuming that all nodes have the 
same chances of being selected) , the probability of all fluxes being different from zero will 
be very high. In this section we wonder whether it happens for all J. 

We will answer this question by computing the fraction of states in the shell having a 
fixed number j of non-vanishing components. If all states are equiprobable, the quotient 
between this number and the total number of nodes in the shell will yield the probability 
distribution of the values j taking into account only abundances of states. 
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We can expect this probability distribution to have a peak at certain value j*. This j* 
will be taken as the typical number of non-vanishing fluxes of the states in the shell. For 
small values of J we know that j* = J. We will see that this is not true for sufficiently 
high J. 

We will now outline the calculation and give the results. The details can be found in 
appendix 1X1 

For any state A G having exactly j non- vanishing components, we define a = j. 
When A is selected at random from with uniform probability, a becomes a discrete 
random variable taking values in the [0, 1] interval whose probability distribution is given 



where Me{j) is the number of nodes in the shell having exactly j non- vanishing compo- 
nents. The formula ( 1571) takes into account only the abundances of states in the shell, and 
hence we are assuming that all states in are equally probable. 

Computation of the quantity NeU) can be achieved using the principle of inclusion- 
exclusion. For simplicity, we will assume equal charges Qi = ■ ■ ■ = qj = q- In the general 
expression fll02p we substitute the number of nodes in the shell and the exact den- 
sity of states ( 1271) . obtaining ( I107p . After normalization, it results in the following exact 
representation for the probability distribution (see ( llOSp in appendix Rl): 



With the assumptions made, we find that -P(a) depends on the radius of the A = sphere 
but it is independent of e. The same method can be used for analyzing the distribution 
P(a) over the whole Landscape, that is, inside the sphere of radius Ri. The resulting 
expression and the subsequent analysis are analogous, and the result is quite similar; in 
appendix [B] the calculation is carried out in the BP regime. 

Using the saddle-point method again, we can approximate the exact formula (158|) by 
(see (I115P in appendix Rl) 

P(a) oc e-^'^") with s{a) = -a log a - {1 - a) log{l - a) + j(f){v, a) , (59) 

where v = q'^s*, and s* is the stationary point of the function 0(s, a) defined in (I58p . The 
saddle point i; is a function of a single variable ha with h = and it has two well-defined 
asymptotic regimes and a crossing regime which requires numerical computation. This is 
plotted in appendix lAj figure |9l 

The distribution ( l59|) has a pronounced peak located at a*{h). This is the typical 
number of non- vanishing fiuxes in the shell (and essentially also in the whole Landscape). 
Its computation must be done numerically by solving the following equation (see (I117P in 
appendix [A]) : 



by 



(57) 



Me 





(60) 
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For each positive value of /i, ( 160|) has a unique solution a*{h) with its own regimes, which 
is plotted in figure HI 



Typical number of non-vanishing fluxes 



o 




1 2 3 4 5 




Sampling the typical number of non-vanishing fluxes 




0.0 0.5 1.0 1.5 




Figure 4: (Left) Numerical solution a*{h) along with its asymptotic regimes shown in equa- 
tions (|119p (small h) and (|12ip (large h). The dashed line is the curve a* = j^. (Right) 
Samples of the typical number of non-vanishing fluxes. The two sampling methods described 
in the text have been used: The inside-shell, maximum-frequency method (blue hollow cir- 
cles) and the inside-ball, average-frequency method (green bullets). The saddle point solution 
displayed in the left panel is also shown (red line). 



Thus, -P(a) is locally Gaussian around its peak 

1 
2 



logP(a) ^ logP(a*) - l-J\s"{a*)\{a - a* f , (61) 



with standard deviation 



a*(l — a*) 1 
, - . , ^ , ^ < . (62) 

Therefore, for large J the peak at a* is very narrow. We can conclude that an overwhelming 
fraction of states in the shell (and in the whole Landscape) have Ja* non-vanishing fiuxes, 
and that for high dimensions this typical number is far from J, which is the typical value 
for the low-dimensional case considered at the beginning of this subsection. 

We should emphasize that the calculation outlined here uses the saddle-point approxi- 
mation, and therefore it is not valid for small J. 

Numerical searches have been carried out varying J for constant q and R, estimating 
a*{h) by counting states. Results are shown in figure H] versus the saddle point curve 
described above. Two sampling methods have been used: 
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• The inside-shell, maximum-frequency method samples states inside a shell and com- 
putes the typical number of non- vanishing fluxes as the value of maximum frequency. 
Some advantages can be mentioned, such as the possibility of performing a better 
sampling of the true set we are describing, but also some disadvantages: the size of 
the sample is smaller, there is an unavoidable intrusion of exterior secant states in 
the shell [231 121], and the sample shows a stronger dependence with the details of 
the lattice. 

• The inside-ball, average-frequency method samples states inside a ball and com- 
putes the typical number of non-vanishing fluxes using the average frequency. When 
compared with the previous method, this method has the disadvantage of sampling 
mostly secant states, but it has the advantage of accounting for bigger sample sizes. 
It also averages over the lattice details, giving a satisfactory agreement with the 
saddle point computation. 

Data used in both samples are a shell of radius R = \/2 and width | with q = 0.15, and 
dimensions J between 2 and 200 in the first method and between 2 and 275 in the second. 
In both cases, the h parameter has been computed using averaged distances, and the lattice 
details can be observed in the first sample in the form of a jagged curve. Note that for 
dimensions up to 7 the typical number of non-vanishing fluxes is J (that is, a* = 1) but 
this abruptly changes to fit the saddle point curve. 

We will close this section with two remarks. First, as the different sampling methods 
considered above suggest, a*{h) curve is very robust as a property of the BP Landscape, in 
the sense that generic subsets of the lattice have a*{h) curves which differ only in subleading 
terms. We can see an example of this feature in appendix [Bl And finally, the role played 
by the exact formula fl271) for ujj{r) in the computation of P{<y) is essential: replacing it 
with the BP estimate (the "pure BP regime") results in a probability distribution valid 
only for h < ^ with the same typical fraction a*{h). Details of this calculation are given 
in appendix [Bl 

3.3 Estimating the minimum positive cosmological constant 

In this subsection we will estimate the explicit dependence of the minimum positive cosmo- 
logical constant with respect to the parameters of the Landscape. We will call A* the actual 
minimum value, and A^ the corresponding estimator. We will assume that all charges are 
equal, for simplicity. In this case, we have 



and we should choose the smallest integer satisfying two conditions: it should yield 
A* > 0, and it should be representable as a sum of J integer squares. If we call such a 




(63) 



i=l 



N 
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number Nj{q), we have the exact formula 



A* = Ao + ^iV,(g). 



(64) 



The computation of Nj{q) can be avoided as we change it by another integer satisfying 
only the first condition: 



Ao + 



> but An 



q\N 



< 0, 



that is, 



2|Ao| 



<N < 



2|Ao| 



+ 1. 



(65) 



(66) 



Thus, is the smallest integer greater than or equal to also called ceiling: 



N 



2|Ao 

g2 



(67) 



If is decomposable as the sum of J squares, we have Nj{q) = N, but in general the 
inequality 

^2|Aor 



< Nj{q) 



(68) 



holds. So we have a lower bound for the minimum value of the cosmological constant, and 
this is our estimator: 



A* > A, 



Ao + ^ 



'^2|Ao 



(69) 



It should be noted that this lower bound is independent of J, even though we can expect 
it to work better when the error for the replacement is small, that is, for large J. In figure 
[5] we show the lower bound along with the actual, brute-force computed minimum, and 
we find a very good agreement for J = 4 and greater. In the figure we can see a straight 
upper envelope of the lower bound, which can be obtained by noting that [x] < 1 + x, and 

thusAo + 4[^l<4. 

Generalizing the preceding argument is not easy, because for different charges we have 

^ = ^o + ^E4^'' (70) 



where q is some average value of the charges. The sum in (1701) is not an integer, and its 
possible values near A = depend strongly on the particular charge values. 

We also lack a reasonable upper bound for the minimum value, because the replacement 
of Nj(q) by Nj{q) + 1 in fl65|) gives a window of width y. The corresponding distance in 
flux space is of the order of the cell spacing, which is too big to be a good upper bound. 
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Minimum A estimate vs. brute-force search (J = 2,3,4) 




I I I I I I r 

-2.0 -1.8 -1.6 -1.4 -1.2 -1.0 -0.8 

log,o(q/R) 



Figure 5: Actual brute-force computation of the minimum positive cosmological constant for 
equal charges in J = 2, 3,4 (hollow circles), along with the lower bound given in (j69p versus 
charge values. The agreement is almost complete for J = 4, and therefore for greater J as 
well. 
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We will adopt another approach from now on. Our strategy will consist of computing 
the number of states of cosmological constant between and and equate it to the 
degeneracy of the minimum positive A states computed in a different way. 

Let A* be a minimum cosmological constant state, that is, A(A*) = A*. There are 
other valid minima of equal value A*. Some of them can be derived from A* using lattice 
symmetries, but there will be a set of minima which cannot be related by any symmetries. 
We will call the set of minima modulo lattice symmetries S*, and its cardinality the 
essential degeneracy. For small J, S* can contain only one state, but we can expect S* to 
grow with J. 

Consider a state A* = ■ ■ ■ ,n}) in E*, and let /fc(A*) be the frequency of the non- 
negative number k in the sequence ■ ■ ■ , |n}|. Note that YlT=o fk — i number 
of different nodes in C that can be derived from A* using the lattice symmetries is 

J! 2-^-^« 



(71) 



The rationale behind equation f l7T|) is to count permutations of the J components except for 
the components which have the same absolute value, and multiply them by the number of 
different "J-quadrants" which can contain such a state (different signs of its components), 
which depend on the number J — /o of non- vanishing components. 

Of course, knowing the sequence {/fc}fc=o,--- ,oo is equivalent to knowing the exact state 
A*, which is a very difficult problem for large J. So we can bound the true value of the 
degeneracy (17T1) between the two extreme cases of states having different values in all non- 
vanishing components (and then, all values are non-degenerate except 0, /fc = 1 if A; 7^ 0) 
and states having all non-vanishing components equal (and then, all frequencies vanish 
/fc = except for two values k G {0, f\ for some €). Then, taking the same number of null 
components /o, 

T"! oJ-fo J] oJ-fo 71 oJ-fo 

/o!(j-/o)! - nr=oA! - /o! ■ ^ ^ 

The degeneracy of the minimum A* is obtained by adding all degeneracies of the different 
states A* G S*, and thus we obtain the bound: 

/o(A*)! (J - /o(A*))! - nr=o /^(A*)! - /o(A*)! ' ^ 

Equating the degeneracy (middle term of f l73|) ) to the number of states in the shell ( 155|) . 
we have 

,4-/o(A*)!(J-/o(A*))! - '-coj{R)J^^, /o(A*)! " ^ ^ 

Note that the shell's width is taken as A^ in flTll) instead of A*, because the number of states 
in the shell is an estimate. Therefore, ([71]) can be taken as a definition of the minimum 
estimator A^. 
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So we are faced with estimating /o(A*), or equivalently, the fraction of non- vanishing 
J 



components a (A*) = \ in terms of which 



^X^) ,4^. [^«(A*)]! - «(A*))]! - - ujj{R) ^4^, [J(l - a(A*))]! 

The sum extended to S* can be replaced by averaging over a probabihty measure of a 
restricted to S*: 

The conditional distribution P(a|S*) is not known, because the set S* is very difficult 
to enumerate. We must assume the robustness of the distribution P{a) computed in the 
previous subsection, and approximate P(a|E*) ^ P{a). We then use the Gaussian nature 
of -P(a) for large J, and we perform the average simply by taking the most probable value: 

R\T.*\ J! 2-^"* 712-^"* 

This equation is not useful because we lack an estimate for |E*|. Nevertheless, for equal 
charges we can replace with (169|1 and use (1771) to obtain a lower bound for the essential 
degeneracy |S* 



UJj{R) [J(l-a*)]! / , q 



-~R J! 2^- I'^^ + y 



^2 



2 A 



01 



q2 



Thus, for the special case of equal charges, we have estimates for the minimum value of 
the cosmological constant and for its essential degeneracy. Figure [6] compares the estimate 
given in fl75]l with actual brute-force computations for low J. 

As said above, generalizing to the case of distinct charges is difficult. Neverthe- 
less, we can assume that for charges not only distinct but incommensurate, the essential 
degeneracy will be = 1. Furthermore, the symmetry degeneracy is reduced to 2'^", so 
that we have the estimate 

Ae - ^ . (79) 

Note that this formula is equivalent to the number of states in the Weinberg Window, f l55p . 
where A/ww should be replaced by the symmetry degeneracy and solved for the width A^. 
But we can improve this by replacing the degeneracy 2*^" by its mean value computed 
using the Gaussian distribution P{q.) given in formulae (|6T|) and (|62|) : 

Jo Jr A/27r«*(l — a*) 

We can check this estimate with the brute-force data for low J (and thus a* = 1) by 
choosing charges with constant geometric mean. In the BP regime u}j{R) depends only on 
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Essential degeneracy E*l estimate vs. brute-force search (J = 3,4) 



o 
o 




iog,o(q/R) 

Figure 6: Actual brute- force computation of the essential degeneracy of the minimum 
positive cosmological constant states for equal charges in J = 3, 4 (circles) along the lower 
bound given in (|78p (thin lines) versus charge values. The strong oscillations of the lower 
bound are caused by the ceiling function (note logarithmic scale). Thick lines are obtained 
replacing the lower bound of by its upper envelope, and they represent an average asymp- 
totic regime which is better followed in the J = 4 sample. 
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Fluctuation of Amin for constant geometric average charge 



1 



Af estimate (J = 2) 
f interval (J = 2) 

estimate (J = 3) 
T interval (J = 3) 



J = 2, q = 0.0248998 



J = 3, q = 0.05371685 



Sample size: 100 

Figure 7: Each box plot represents a sample of 100 choices of J charges with constant 
geometric average charge for which Amin has been computed by brute-force search. Thick 
lines correspond to estimate (I79p . which does not fluctuate (at least in the BP regime). These 
lines cross inside the boxes in good agreement with numerical data. Dashed lines enclose the 
values provided by the "thermodynamic" estimator T* , which show even better agreement. 
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the geometric mean and therefore it does not fluctuate. Figure [7] shows that (!79|) is a good 
estimator for the mean value of such fluctuations. 

We can derive another estimator with an additional parameter which allows some con- 
trol of the bias by providing an interval instead of a single value. Nevertheless, the interval 
is not a confidence interval and we lack a rigorous proof of strict inclusion. 

We begin by considering the following partition function of the Landscape: 

Z = Y,e[m]e-'^ . (81) 

Here, the sum is carried out over all lattice nodes, but the step function excludes all 
negative cosmological constant states. A plays the role of the energy, and T is the associ- 
ated temperature. Low lying states constitute the dominant contribution to the partition 
function at very low temperatures, that is, 

Z^^e^*e~^+eV^ + --- , (82) 

where A* is the actual minimum positive value of the cosmological constant of degeneracy 
, and A is the first excited state of degeneracy . Introducing the gap (5A = A — A*, 
the low temperature free energy is 

F=-T log Z A* -TS* -T log (l + 6^-^*6"^ + ■■■), (83) 

and, at sufficiently low temperatures, namely T ^ 6A, we have the linear dependence 

F ^A* -TS* . (84) 

At first sight, this equation can be used to estimate the ground state energy A* and its 
entropy S* by computing the free energy in another way. For this, we use the counting 
measure in the Landscape flj{r) given in f l29l) . and write 

Mr) 

Z= e"— dfij(r), (85) 
Jr 

where we have used the step function to cut off the interval of integration. A depends on 
the radial variable r in flux space: 

A{r) = ^-{r' - R') . (86) 

Note that the counting measure is reminiscent of the spherical volume element (divided 
by volQ) in J dimensions, and in fact reduces to it in the BP regime, but it is really a 
discrete distribution whose properties are crucially different from the continuous ones, as 
will be seen below. 
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Replacing df2j(r) = ijjj{r)dr and using the exact contour integral (127|) for uj{r) we 
obtain, upon reversing the order of integration and changing the variable r to m = A(r), 



e T ujj[r) ar = 

R JR 



1 
1 

27Ti 



I -Ji i=i J 

e--^+''-^2r dr \ JJ ^{q^^s) ds (87) 
i=i 

/"OO ^ J 

2 / e-^^'^'^''du\e''^"\\§{qls)ds, 
Jo J 



and the radial integral converges provided Re{s} < that is, if the pole at ^ is located 
to the right of the integration contour 7. Then, we have a contour integral representation 
for the partition function: 



We can see that flSSl) gives the initial expression fIST]) by expanding the {} sums under the 
product sign. Interchanging sum and integral, and evaluating each integral, it turns 

1 r e^R'-sj:Unhl 



■ds. (89) 



If the factor multiplying s in the exponent is negative (thus A is positive), then we can 
close the contour by a large half circle to the right; the integral on the circle vanishes, and 
the contour encloses the pole at with negative orientation, resulting in a residue e~T. 
On the other hand, if A is negative, the contour can be closed on the left, but the integrand 
has no poles in this region, and therefore the integral vanishes. Thus, we recover (IHTI) . 

Equation ( l88l) provides an independent evaluation of the partition function by numerical 
computation of the contour integral. Then we can use the values obtained in the low 
temperature region to fit equation ( |84l) and estimate A* and S* . Nevertheless, the expected 
straight line (|8^ is not obtained this way, but a rather different behavior, as we now explain 
(see figure [8]). 

Assume that we want to approximate expression fl85l) at low temperatures. Changing 
the integration variable to m = A(r) again we obtain 

Z = I e-^ojj{r)dr = I e^^wj f Vi?^ + 2u) ^=^= , (90) 
jR Jo ^ ^vR'^ + 2u 

which expresses Z as the Laplace transform of If this function is continuous 

near u = 0, the asymptotic behavior of Z is simply 

Z^I^I^, (91) 
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Free energy 



Entropy 
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Discrete (linear app.) 
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Numerical 



Discrete (exact) 

— Discrete (iinear app.) 

Continuous 

o Numericai 



Oe+00 



2e-05 



4e-05 



6e-05 



Oe+00 



1 1 1 

2e-05 4e-05 6e-05 
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Figure 8: Comparison between thermodynamic magnitudes preserving discreteness or assum- 
ing continuity of the density of states. Near T = 0, the exact free energy is locahy hnear (left 
panel), and the exact entropy reaches the value log (4) (right panel). Note that the numer- 
ical computation of the contour integral (|88p follows the continuous curves in both panels. 
The exact entropy (solid thick blue line) and the continuously approximated one (solid thin 
red line) never intersect, but the linear approximation (green dashed horizontal line) does, 
defining the estimator T* given in (I95p . 
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and the thermodynamic magnitudes are, when T — )■ 0, 

Tuj{R) 



F = -T log ■ 
^ = 1 + log ■ 



R 

TujjR) (92) 
R ' 
{A) = F + TS = T . 

As seen in figure [8], the numerically computed partition function is a good approximation 
for curves computed assuming continuity, suggesting that the numerical quadrature rule 
works as if the density of states were continuous, and therefore it is not useful to estimate 
the ground state of the discrete spectrum. A* and S*. 

It is worth emphasizing the different behavior of the thermodynamic magnitudes be- 
tween the discrete and continuous systems. Any computation method for Z should respect 
the discrete nature of the spectrum. Otherwise the method would accurately compute the 
resuh (|92|). 

Nevertheless, the two behaviors coincide at high temperatures, suggesting that there is 
a crossing temperature T* below which the difference becomes drastic. Of course, in figure 
M (right) we can see that the entropy of the discrete system (obtained by exact computa- 
tion of Z using brute-force compilation of states) and the approximate one (that of the 
continuous system) never cross. But if we approximate the former using the approximated 
free energy (183!) . then we can find a temperature at which the continuous and discrete 
entropies coincide. This temperature T* signals the point of a drastic deviation of the 
discrete system from the continuous one, and it will be interpreted as an estimator for the 
ground state A*. 

Thus, we need the entropy in ( l92l) . which will be termed ^cont, and the entropy derived 
from (!83l) . which is 

^disc = S*+\og(l + e^-^*e-t) + . (93) 



T(l + e^-^*e-¥) 

The crossing temperature is defined by the equation S^isc = 'S'cont, but fl93|) contains the 
degeneracy entropies S*, S and the gap 6 A as additional parameters. Of course, we cannot 
use a single equation to fix four parameters; we will compute the degeneracies as 2^^"* for 
the case of incommensurate charges (that is, S* = S), and we will introduce the parameter 
= P", so that the equation for T* reads 

S* + log(l + e-") + = 1 + log T^A^ (94) 

and gets solved as 

T* = e'*^-(l + e-^)eTf^-K (95) 

In terms of the previous estimator A^ ( I79|) . we have 

T*{r]) = f{r])K, with /(r/) = (1 + e-") eirW-i , (ge) 



28 



where the prefactor /(?7) is of order 0{1) for < rj < oo: it is a monotonically decreasing 
function with /(O) = 2e~^ and f{oo) = e^^. Thus, when we let rj run across its range, the 
estimator T*{ri) spans the interval [l,2]e^^Ae instead of the single value A^. Nevertheless, 
we cannot give an argument which favours the inclusion of the true value inside this interval. 
Despite this uncertainty, our numerical searches validate the given interval for J = 2 and 
J = 3: it contains the median and its width is smaller than the interquartile range (see 
figure [7]) . 

Thus, we can conclude that T* is a good estimator of the charge-averaged minimum 
cosmological constant. 

3.4 A possible influence on the KKLT mechanism 

The exposition in this paper is restricted to the Bousso-Polchinski Landscape, which is an 
oversimplification of the string theory Landscape. The Giddings-Kachru-Polchinski model 
[26] is a more realistic approach to the true Landscape, and it can be endowed with a 
mechanism for fixing the moduli of the compactification manifold, the so-called KKLT 
mechanism [27]. In this model, the compactification moduli are fixed by the presence of 
fluxes and corrections to the superpotential coming from localized branes. The fixing of 
the moduli lead to metastable dS states by the addition of anti-branes and exceeding flux. 
The model can be further corrected to yield inflation of the noncompact geometry by the 
repulsion between a brane and an antibrane, both located at a Klevanov-Strassler throat 
[28] . The important point is that both moduli fixing and brane inflation need the presence 
of flux quanta. 

As far as we know, there is no combination between the BP Landscape and the KKLT 
mechanism, in the sense that there is no known realistic model in which all moduli are 
fixed and a large amount of three-cycles lead to an anthropic value of the cosmological 
constant without any fine-tuning at all. Nevertheless, there are some toy models as the 
six-dimensional Einstein-Maxwell theory [291 EQ] in which it is possible to identify a Land- 
scape with all moduli fixed, and it has a dual model with fiux; or the more sophisticated 
(yet unrealistic as well) models coming from F-theory fiux compactification [31j or IIB 
compactifications with fiuxes [32], to name a few. Thus, it is plausible that a complete, 
realistic model exhibiting all these features will be built in the near future. 

Let us assume that such a model exists. Furthermore, let us assume that the a*{h) 
curve discussed in the section [3^ can be generalized, in the sense that, as a characteristic 
feature of this conjecturally realistic model of the string theory Landscape, there is a typical 
occupation number of the fiuxes that is generically different from 1 and that vanishes when 
the number of fiuxes J is too large. Then, we should conclude that, generically, there will 
be a finite fraction 1 — a* of three-cycles with vanishing fiux, which may be dominant if J 
is large and the charges qi are not too small. This fraction of vanishing fiuxes can spoil the 
stabilization mechanism of the corresponding moduli and can also spoil the brane infiation 
scenario if the three-cycles located at the tip of the KS throats are devoid of fiux. 

Thus, the a* fraction we have found in the BP Landscape can be an obstacle for the 
commonly accepted KKLT mechanism if it is also present in a realistic Landscape model. 
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4 Conclusions 



We have developed an exact formula for counting states in the Bousso-Polchinski Landscape 
which reduces to the volume-counting one in certain (BP) regime. The formula might be 
useful in BP Landscape examples when the parameter h — ^ is large enough to invalidate 
the BP count. Applications of the formula are given which avoid volume-counting: counting 
low-lying states, computing the typical fraction of non- vanishing fluxes or estimating the 
minimum value of the cosmological constant which can be achieved in the BP Landscape. 
Numeric computations and brute-force searches have been carried out to check the results 
of our analytic approximations, and we have found remarkable agreement in all explored 
regimes. 

In particular, we have discovered a robust property of the BP Landscape, the curve 
a*{h), the typical fraction of non- vanishing fluxes, which reveals the structure of the lattice 
inside a sphere for large J as the union of hypcrplane portions of effective dimension near 
J a*. This result is important in computing degeneracies, which are used in estimating 
the minimum cosmological constant. We have not developed a formula for this minimum, 
not even a probability distribution, but rather an estimator which can produce acceptable 
results in predicting the mean value of the minimum for fluctuating different charges around 
its geometric mean, or essential degeneracies in the case of equal charges. 

Finally, we have pointed out that if we can generalize the typical number of non- 
vanishing fluxes in a realistic model describing the string theory Landscape, then it could 
be an obstacle for the implementation in the same hypotetical model of the KKTL moduli 
stabilization mechanism. 

Thus, the exact formula improves the way of counting as compared to previous propos- 
als, and should be considered in all those problems which require state enumeration in the 
BP Landscape. 
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A Detailed computation of the typical fraction of non- 
vanishing fluxes 

Consider an arbitrary subset of nodes in the lattice S C £. We will decompose its cardi- 
nality I S I = As as 

J 

Afj: = J2^^^3), (97) 

j=0 

where MsH) is the number of states in S with exactly j non- vanishing components. If a 
is the fraction j then its probability distribution over S is 

This formula takes into account only abundances of states in S, and hence it assumes that 
all states in E are equally probable. 

Let us introduce some notation. We will call the set of indexes of components J = 
{1, 2, ■ ■ ■ , J} and L, M will denote any subset of JT". The symbol T,l will denote the set 
of states of S having (at least) vanishing components outside L or, in other words, the 
intersection between S and the subspace spanned by the directions in L. Ml will be the 
number of elements of S^. Thus, T^j comprises all E states, that is, Nj = Afs, and S0 
only contains the node at the origin, so that A/0 = 1. The inclusion-exclusion principle in 
Combinatorics allows us to write 

•^^(^■) = EE(-i)'"' E-^^- (99) 

Lcj 1=1 MCL 

\L\=j \M\=e 

The idea behind ( 199|) is that we cannot simply sum all states lying inside every E/, G S, 
that is, equation A/E(j) = Y]lcj -^l is not true, because of the intersections between 

\L\=3 

different subsets S^. The nodes inside these intersections are removed twice, and they 
should be added again, and so on. This is the cause of the alternating sign in f l99|) . 
Now, we must deal with expressions of the form 

E E -^^ ■ (100) 

LcJ MCL 

\L\=j \M\=e 

All the subsets M d J occur in the previous sum the same number of times, which 
coincides with the number of supersets L inside J (with \L\ = j) of a given M (with 
\M\ = i < j). This can be computed as the number of subsets of J\M of exactly j — I 
elements, that is, 



ee^«=(^:')5:a^«- 

LdJ MCL ^ MCJ 

\L\=j \M\=l \M\=i 



(101) 
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By substituting (fTOB in ([99]), we get 



(102) 



\M\=i 

We can check that equation (11021) satisfies the normahzation condition ([9 

J j 



i=o j=o e=i ^-^ ^ Mdj 



E 

£=0 



M 

J 

\M\-- 



MCJ 



(103) 



V . ^\M\-- 

MCJ 
\M\=J 

Now, we will take the S set as the states of the lattice C near the A = surface, inside a 
thin shell of width A^, whose number is (llOOp 

Ar, = ^^A.. (104) 

With different charges gj, different subsets M G JT" with the same cardinality \M\ = i will 
not have the same number of states N'm- However, in the simplest case where all charges 
are equal, all the number of states coincide: 

E-^"=(')^A., (105) 



MCJ 
\M\=i 



where we have used that there is ('^) different subsets M G JT" with \M\ = £. Substituting 
(llOSp in (I102p and reordering the binomial coefficients we have 



R\J. 



(106) 



^.o)^E(-ir'(;:;)(^)^A. 
:)e (9 

Now, we can substitute the exact integral representation specialized for equal charges 
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for uJe{R) and perform the binomial sum: 



AefJ\2R 



r J 



ds 



2A, // 



(107) 



27ri \j 

Normalization provides the probability distribution of a = j we are looking for: 

'2R f J\ I 



) — [ e'^(''") ds with 0(s, a) = sR^ + aJ\og[d{q^s) - l] . (108) 



The next step is estimating P{a) using the steepest descent method again. The equation 
for the saddle point is 

2 ^\<fs) 



b\s) = R^ + aJq^ 



0. 



(109) 



d{q'^s) - 1 

As we did before, we can find approximate expressions for the saddle point in the two 
regimes of d function. If s* is the saddle point, convenient variables are v = q^s* and 
h — in terms of which we have 



R^ r 

0(v) ^ ^ v + ha\og[e(v) - 1] 
In the large v regime, we have -^Iv) ^ 1 + 2e~^ + 2e~^'" ^ so that 

+ ha{-v + log2 + log(l + e-2")) 



(110) 



R^ 






V 










V 





(111) 



and the corresponding saddle point equation has a solution 
In the small v regime, we simply have 'd{v) ~ y^, so that 



(112) 



q^ 
R^ 

n1 



V + ha 

V + ha 



I log TT - i log + log ( 1 - 
|log7r - ilogt; - 



V _v_ 

TT 27r 



(113) 



The corresponding saddle point equation is quadratic in -j^^, and its solution is 

47r 



v{ha) 



ha 2 {ha)i 3{ha)^ 



(114) 
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Saddle point of the probability distribution 




Both approximate solutions (11121) and (11141) . as well as a numerical one, are plotted in 
figure m Note the logarithmic divergence at ha = 1 of the saddle point. For larger values, 
v{ha) becomes complex, and the integral P{a) begins to oscillate, losing its probabilistic 
meaning. 

Using the asymptotic form of the binomial coefficient, we can write 

P(a) oc e-^'^") with s(a) = -a log a - (1 - a) log(l - a) + ^0(t;, a) , (115) 

where (j){v, a) also depends implicitly on a through v. The exponent s{a) has a maximum 
a*, given by the equation 

ds d(j) dv d(f) 

— = - log a + log(l - a) + — -J- + ^ 
da .ov, da oa 



;ii6) 



1 — a 

log— — + \og['&{v (ha)) - 1] = 



In the first equality of (I116P we have used the definition of the saddle point v, and in the 
second we have used equation (IllOp . We obtain a* as the unique real solution in the [0, 1] 
interval of 

^\v(ha)] = ^—. (117) 
-' 1 — a 

Note that the saddle point v{ha) is defined only for ha < 1 (see fig. [H]), so that the 
right hand side of equation (I117p as a function of a has domain a G [0, In this interval, 
^[viha)] decreases from oo to 1, while the right hand side increases from 1 to oo in a G [0, 1]. 
These conditions guarantee the existence of a unique real solution a*{h) G [0,1] for all 
positive h. Figure [TO] shows the construction of a*{h). 

For small /i, we can use the small ha regime for v{ha) (eq. (I114p ) and -i? (eq. ( I32|) ). in 
which equation (11170 reads 

ha 1 — a 
and has the solution 

a*{h)^l-—\^/h(hT8^-h]. (119) 

4:71 

For ha approaching 1, we can use (11120 for v{ha) and the corresponding regime for ^, so 
that (I117P simplifies to 

l + 2e^'°^K7iTJ-i) = (120) 

1 — a 

We can isolate h{a*) in (11200 . obtaining 

h{a*) ^ ^ , (121) 

which for large h is simply a* ~ -i. All these regimes are shown in figure H] along with the 
numerical solution of equation (I117p . 
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Construction of a (h) (0 < h < 3) 




Figure 10: Construction of the solution a* of equation (jll7p for < h < 3. Dashed curve 
corresponds to h = 1. For h <^ 1, the solution is near a = 1, and for h ^ 1 the solution 
approaches zero. 
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B Typical fraction of non- vanishing fluxes in the pure 
BP regime 

We can repeat the computation of the probabihty distribution of the fraction of non- 
vanishing fluxes -P(a) (given in appendix |A]) in the pure BP regime step by step. Using 
the previous notation, we will consider all states inside a sphere of radius r as the set 
S C £, but we will use the BP formula ([7]) instead of the exact integral representation of 
the number of such states (hence the name "pure BP"): 

U^ = ^j{r)^ (122) 

jr(f) voiQ 

The computation of the probability distribution is based in the decomposition (p7|) . which 
we repeat here for convenience: 

J 

Af^ = J2^dj), (123) 

j=Q 

After using the inclusion-exclusion principle, we have f ll02p 

^f^iJ) = E(-ir'(^:n E (124) 

i=l ^ MCJ 

\M\=i 

As before, we will assume that all charges are equal vol Q = . Thus we have an equation 
analogous to f llOSp 



\M\=t 

where we have used the parameter 



i = (126) 



Plugging fll25p in the decomposition fll24p . 



Now we take the Hankel definition of the gamma function: 

1 1 



Tix) 2Tii 



e'z-'^dz, (128) 
c 
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where the contour C in complex z plane encloses the origin coming from — oo below the 
negative real axis, turning around the origin and leaving towards —oo over the negative 
real axis. By substituting 



r(| + l) 2m J c 



(129) 



in ( I127p . we can interchange the sum and the integral and the binomial theorem allows us 
to write 



jj2TTiJc \^/z J Z 

By changing the z variable to s = zj^ the contour only gets scaled and can be deformed 
back to its original form, resulting in 



We have defined the function 



c 



0(s) = ^ + «log(^-^-lj, (132) 

which depends on the two parameters mainly used in the remainder of the paper: 

J cj^ j 
h = — — and a = — . (133) 

J 

Note that ^ = ^J^. 

We can normalize A/e(j) dividing by A/s in f ll22p . thus obtaining the following integral 
representation for the probability distribution: 



(134) 



We can evaluate the integral using the steepest descent method, assuming that J is large 
and that h and a are constant parameters. Furthermore, h cannot be too large for the 
saddle point approximation to remain valid. 
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Thus, for large J, -P(a) is 

P{a) oc e-^'^'^'^ with s{a) = -a log a - {1 - a) log(l - a) + 0(s*) , (135) 

where we have used the asymptotic expression for large J of ( and s* is the saddle point 
of 0(s) which is compatible with the integration contour C. We compute s* as a solution 
of 

, , , 7r 1 , , 



which can be rewritten as 



Sy/s-S + P = 0, (137) 



with the new parameter = |f • Let us assume for simplicity that u is small (that is, h is 
small); then f ll37p has two solutions in the positive real axis, one of them near and the 
other near 1. Near s ^ 0, 4>"{s) ^ ^ > and near s ^ 1, (i)"{s) ""^^^ < 0, so 

the first corresponds to a local minimum over the positive real axis and the second to a 
local maximum. But the integration contour crosses the positive real axis vertically from 
lower to upper half plane, and therefore the integrand has a maximum along C if it has a 
minimum along the real axis. So our saddle point is near for small and it follows from 
(11371) that we must have 

s*^u. (138) 



We can find the maximum of the "entropy" s{q) in (]135p . which we will call a*, by solving 
the stationary condition 



ds(a) _ ^ g ~ " ' ^^^^^ '^^* 1 ^^(^) 



da a ds 



=s*, da da 



log^^ + logf^-1 ) =0. (139) 



s=s* a 



With the estimate fll38p for small eq. (I139p can be rewritten as 



a 



a 



V. (140) 



In the preceding equation, if v is small, then its solution a* must be near 1; replacing a ~ 1 
in the denominator of (I140p we obtain exactly (IllSp . which determines a*(/i) for small h. 
So this approach ends up with the same estimate for the typical fraction of non-vanishing 
fluxes for small /i, despite the fact that we are computing this fraction a* over the whole 
Landscape instead of over a thin shell. This is a strong indication that the result obtained 
for a* is robust, in the sense that it does not change significantly for different generic 
subsets of the Landscape. 

We remark that the robustness of the a*(/i) curve is not a feature of spherical symmetry. 
This is exemplified by the set of secant states, which is not spherically symmetric [23| [2l]. 
If we repeat the computation in this appendix taking the set of secant states at distance 
r as the subset S C £ instead of the ball of radius r, we would obtain the same estimate. 
Thus, we should consider a*{K) a robust property of the BP Landscape. 
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The preceding computation used the BP regime (11221) . and we have mentioned above 
that h should be small for the validity of the saddle-point approximation. Now we can 
wonder, how big can h be before invalidating the approximation. Can we somehow continue 
this result to higher values of hi 

The answer to this latter question is negative, as can be seen in figure ITTl which, in turn, 
answers also the former question. In figure [11] we show the saddle points of as they 



Evolution of the critical points of ^(s) witfi v 
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Figure 11: Evolution in complex s-plane of the critical points of the cf) function (jl32p . exponent 
of the integral representation of P{a) in the pure BP regime, versus the v parameter. The 
saddle point s* plotted in red (which begins at s = and enters the lower half plane for 
V > ^) is compatible with the integration contour C, which means that C can be deformed 
to coincide with the steepest descent contour crossing trough s*. 



evolve in complex s-plane from = to z/ = 4. Both of them start at s = 0, 1 respectively 
for z/ = and evolve along the positive real axis until they coincide at s = | for z/ = ^. 
This coincidence can be seen by writing (11371) in a variable x = y/s and considering its 
derivative, which is a polynomial whose roots x* are the critical points of p{x): 



p{ 



x] 



Q 2 

X —X +Z/ 



d 

dx 



3x-2x = 



X e 



{o,^} ^ P{xle{u,u-^], (141) 



that is, p{x) has two critical points with values above and below the x-axis only ii u < ^. 
For u > both critical points have positive values, and the roots of p{x) leave the real 
axis. 
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Only the root near is compatible with the integration contour C, that is, C can be 
deformed to meet the steepest descent contour crossing by the saddle point near s = 0. For 
1/ > ^ both roots leave the real axis; as a consequence, the integral -P(a) remains real but 
becomes non-positive and rapidly oscillating, thereby losing its meaning as a probability 
distribution. Therefore, the pure BP regime is only valid for u < that is, for 

/i< — = 0.93084. (142) 
27 ^ ' 

This value represents the upper limit of validity of the pure BP regime in the computation 
of P{a). We can intuitively understand the existence of this limit by considering that the 
inclusion-exclusion principle f ll24p . being an alternating sum, is very sensitive to inaccurate 
computations of the cardinals of the subsets appearing in the sum. Therefore, when we 
use the BP formula (11221) for computing the cardinals of the subsets with growing h, that 
is, for decreasing r, the error in the formula causes strong oscillations in P{a). 
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